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Abstract. 

In the hard pomeron model consequences are studied which follow from the recently 
obtained form of the diffractive amplitude for the double scattering on the nucleus and the 
related EMC effect at small x. It is shown that at large Q 2 the double scattering contribution 
to the latter falls as Q -0 - 6338 and in all probability dominates the total effect. 
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1 Introduction. Double scattering 

Perturbative QCD predicts that hadronic scattering at high energies proceeds via the ex- 
change of a hard (BFKL) pomeron, with an intercept above unity. Much effort has been 
spent in analysing possibilities to see a signature of the hard pomeron in the observed be- 
haviour of the structure functions at low x and in jet production. 

Recently a new object has drawn considerable attention within the perturbative approach, 
namely, the diffractive amplitude and the triple pomeron interaction related to it [1-4]. The 
found diffractive amplitude turned out to have a rather specific form. Since this amplitude 
governs the double scattering contribution to the cross-section on the nucleus, these find- 
ings open new possibilities to compare predictions of the perturbative approach with the 
experiment. 

In short, the diffractive amplitude calculated perturbatively and in the high-colour limit 
has been found to possess the following properties [1,4] 

1. The double pomeron exchange contribution to it can be completely absorbed into the 
triple pomeron one. 

2. At small momentum transfers k the latter is singular. 

3. At small k the coupling of the pomeron to a hadronic projectile behaves as Itlk 
These properties have been established in a pure perturbative approach, with a fixed 

strong coupling constant a s . It is not known in what manner they will change if one introduces 
confinement effects in any form. 

In this note, rather than speculate on the latter problem, we study the observable con- 
sequences of these properties for the double scattering on the nucleus. To stay within the 
perturbative appproach the projectile has to be either a heavy "onium" or a highly virtual 
photon. Having in mind applications to the low x EMC effect, we choose the latter choice. 

The amplitude A2 corresponding to the double scattering on the nucleus and normalized 
according to 2 ImA = a tot is given by the well-known expression (for A > 2) 



where pa is a two-nucleon density. If it is taken factorized (no correlation approximation) 
then 




(1) 



Here Fa(k) is a two-nucleon form-factor of the nucleus: 




(2) 



F a (k)=TI(k) 



(3) 
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where Ta(k) is a two-dimensional Fourier transform of the standard nuclear profile function 
7a(£>). For the deuteron, instead of (2), 

F 2 (k) = J d 2 r\^(r)\ 2 e™\ k z = (4) 

where ip is the deuteron wave function. 

The high-energy part a 2 (k) is given by an integral 

1 f°° 

ia 2 (n) = ^2 J Q ds 1 Disc Sl D(s,s 1 ,n) (5) 

where D(s,si,k) is a diffractive amplitude for the cm. energy squared s, diffractive mass 
squared s\ and transferred momentum k. 

We recall that the eikonal (Glauber) approximation consists in taking in (5) only a con- 
tribution from the intermediate state equal to the initial projectile hadron, which gives a 2 
independent of k: 

cl 2 (k) = a 2 (6) 

where a is the forward projectile-nucleon amplitude (also normalized to 2Ima = <j tot ). Then 
(1) and (3) immediately lead to the well-known result 

A 2 = iC\a 2 w 2 , w 2 = J d 2 bT 2 A {b) (7) 
2 Triple pomeron contribution 

The diffractive amplitude found in the perurbative QCD is totally given by a contribution 
from the triple pomeron interaction. So the corresponding double scattering amplitude on 
the nucleus will have an essentially non-Glauber form. 
The discontinuity entering (5) turns out to be given by 

Dl ee iDisc^M^) = ^N 2 (N 2 - 1)£ jf[d\ r -^ 

exp(-iK(r 2 + r 3 )/2)(p 1 (s 1 ,0,r 1 )(p 2 (s 2 ,K,r 2 )(p 2 (s 2 ,-K,-r 3 )5 2 (r 1 +r 2 + r 3 ) (8) 

Here 4>i(s, k, r), describe the upper (i = 1) and two lower pomerons (i = 2, see Figure) for the 
energetic variable s, momentum k and the intergluon transverse distance r. N is the number 
of colours; s 2 = s/s\. The pomerons are joined by the triple pomeron vertex, whose form 
was obtained in [4] in the high-colour limit. Presenting the <5-function as an integral over the 
auxiliary momentum q we rewrite (8) in the form 
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where 



and 



Xi(s,0,q) = J d 2 rr 2 V 4 (j)i(s,0,r)expiqr (10) 

(ii) 

The solutions 1 ( 2 ) can be obtained by using the Green function of the BFKL equation 
for a given total momentum Gi(s,r,r'). For the projectile (see [5]): 

r) = | dV(G s ( K ,r,0)-G s ( K ,r,r') Pl (r') (12) 

Here pi(r) is the colour density of the projectile as a function of the intergluon distance with 
the colour factor (l/2)5 a b and g 2 separated. As mentioned, for the projectile we choose a 
highly virtual photon with p\ = —Q 2 < 0, which splits into qq pairs of different flavours. The 
explicit form of p is well-known for this case [5]. For the transverse photon 

2 N f 1 

P?V) = ^E 2 // da(m 2 K 2 (e f r) + (a 2 + (1 - a) 2 )e 2 K 2 (e f r)) (13) 
j J 

where = Q 2 a(l — a) + m 2 and m/ and Zf are the mass and charge of the quark of flavour 
/. For the longitudinal photon 

2 N f i 

p[ L \r) = ^Q 2 J2zj [ daa 2 (l-a) 2 K 2 (e f r) (14) 



For the hadronic target, we assume an expression similar to (12) with a colour density P2{r) 
non-perturbative and its explicit form unknown. For our purpose it is sufficient to know that 
the corresponding mass scale is not large. 

The two lower pomerons are in their asymptotic regime. So \2 can be found using an 
asymptotic expression for the pomeron Green function, in the same way as in [6], to which 
paper we refer for the details. One obtains then 

X2(s,K,q) = 8s A (it /p In s) 3 / 2 F 2 (k) J ( K ,q) (15) 

where A = (a s iV/7r)4m 2 is the pomeron intercept, (3 = (a s iV/7r) 14^(3), 

JU,q)= I — —, r-4 n r ( 16 ) 

V ,HJ J 2ir \k/2+p\\k/2 -p\\q+p\ 

and F2(k) desribes the coupling of the lower pomeron to the target 

dRdr \R + r/2\\R-r/2\ = 
At small k F 2 diverges logarithmically: 

F 2 ~ -4(R/N)]nn, n^O (18) 
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where 



N > A 



R = — I d z rrp 2 (r) (19) 



has a meaning of the average target dimension. 

If one also takes an asymptotic expression for the upper pomeron then the whole discon- 
tinuity D\ aquires a factorized form [6] corresponding to the triple pomeron contribution in 
the old Regge-Gribov theory: 

£>i(a,ai,/c) = --f 1 -fi(K)j 3P (K)P(s 1 ,0)P 2 (s/s 1 , K ) (20) 

Here P(s, k) is the pomeron propagator 

P(s, k) = 2^s 1+A {(3 In s)~ e (21) 

with e = 1/2(3/2) for k = 0(> 0). The vertices 7 1 / 2 ) describe the pomeron interaction with 
the projectile (target). Using (13) and (14) one obtains, neglecting the quark masses 

ln 2 Q s 



7i = a s VW^T(b/Q)exp (22) 

where Q 2 is the photon virtuality and for a transverse (T) and longitudinal (L) photons 

9vre 2 Z 2 , 2, 

b T = , b L = -br 23 

256 9 

Z 2 = J2f Zf- The target vertex j 2 is proportional to F 2 , Eq. (17) 



72 (k) = a s VN 2 - 1F 2 (k) (24) 

and so logarithmically divergent at small k. The triple pomeron coupling 73P is also singular 
at small k: 

32a 2 N 2 B 

wiK) = TWTt^ (25) 

where the number B is given by a 6-dimensional integral over auxiliary momenta [4]. Nu- 
merical calculations give 

B = 4.98 ±0.01 (26) 

in agreement with the analytic result recently obtained in [7,8] (their g^p = (2tt) 4 B). 

Putting explicit expressions for P's and 7's into (20) and performing the integration over 
si one finds 

R 2 s 2A ln 2 K 

a2 = ~ C2 ^ c^f^t (27) 

where c 2 is a known number. From this purely asymptotic expression one could come to 
two conclusions. First, the double scattering contribution behaves in Q exactly as the single 
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scattering term, that is, as 1/Q. As a result, the EMC effect related to it should be indepen- 
dent of Q. Of course, this is an immediate consequence of the factorization property of the 
asymptotic triple pomeron interaction. Second, the singularity at small k leads to a stronger 
dependence on A. Since roughly speaking k ~ 1/Ra where Ra is the nuclear radius, from 
(27) one could conclude that the double scattering is enhanced by a factor (Ra/R) In 2 (Ra/R) 
as compared to the standard eikonal result. 

However all these conclusions are in fact wrong, since the upper pomeron enters the triple 
pomeron interaction not at asymptotic energies a s lns » 1 but at lower ones a s lns ~ 1, 
where the asymptotic expression for its Green function used in (20) is not valid. So we have 
to recur to the exact expression for it. Due to azimuthal symmetry of the projectile colour 
density we can retain only terms with zero orbital momentum in it: 

■co (^ 2 + 1/4) 2 {r/r ' y2W ^ 

where 

u(v) = 2(a s iV/27r)(Y>(l) - Re^(l/2 + iv)) (29) 

If we take this expression, put it into %i, Eq. (10) and integrate first over s±, as indicated in 
(10), and afterwards over r, we obtain 

dsisr 1_2A x(si, 0, ?) = -—/ d 2 r'r' Pl (r')I(q, r') (30) 



where I(q, r') is the remaining integral over u: 

vy ' 1 J 2A-lo(v) r(l/2 + ii/ v ' 

The asymptotic expression discussed above follows if one takes in (31) all terms except 
the denominator 2 A — uj(v) out of the integral at v = and in the denominator presents 
uj(u) = A—/3i> 2 . Evidently this procedure is wrong. In fact, the integral (31) can be calculated 
as a sum of residues of the integrand at points v = ±ixk, < x\ < X2 < at which 

2A-w(^) = (32) 

Residues in the upper semiplane are to be taken if qr'/2 > 1 and those in the lower semiplane 
if qr'/2 < 1. Thus we obtain 

^yn^E^V/^ (33) 

where 



(±) T(l/2Tx k )/T(l/2±x k ) 

k ^(l/2-x fc )-V'(l/2 + x fc ) 1 J 



and the signs should be chosen to always have {qr' /2) ±2xk < 1. 
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The first three roots of Eq. (32) are 

si = 0.3169, x 2 = 1.3718, x 3 = 2.3867 (35) 
with the corresponding coefficients 

c{ +) = 0.1522, 4 +) = -0.1407, 4 +) = 0.03433, 

c {-) = 0.007866, 4 _) = -0.001802, 4 _) = 0.004494, (36) 

Returning to Eq. (5) for a 2 as an integral of the discontinuity D\ and putting expressions 
for X2 and the integrated X2, Eqs. (15) and (30), into the latter we obtain 

a 2 = -128^iV(iV 2 - l)^_^|( K )]r J d 2 rr Pl (r) 



where 



i#\l) 2xk Bl + \r, K ) +cP {-)-** Bi~\r, K )] (37) 

^ ±)(r ' K) = / 7^^ fcj2 (^9-K/2)e(±(2/r-g)) (38) 

As we observe, in the general case the factorization property is lost: the integrals B^' 
depend nontrivially both on the projectile and target variables. However one can see that 
this property is restored in the limit of high Q 2 , relevant for the hadronic structure functions. 
In fact, in this limit the characteristic values of r are small: r ~ 1/Q- Let us study how 
the integrals B^ behave at small r. In B^ evidently large values of q are essential. The 
integrals J behave as Inq/q at q — > oo. This leads to the following behaviour. 

(r, k) ~ r 1_2:rfc if 2x fc > 1 and ~ const if 2x k < 1 

s£ _) (r, k) ~r l+2xh 

Combining this with other factors depending on r we see that all terms multiplying p\ in the 
integrand behave as r 2 at small r, except the first term with k = 1, which, due to 2x\ < 1, 
behaves as r 1+2xi . Evidently this term gives the dominant contribution in the limit Q 2 — > oo, 
when (37) simplifies to 

a 2 = -^Va^iV 2 - l)^ J ^^n) (39) 

where the numbers b\ and B\ are very similar to our old b and B with additional powers of 
the variable in the integrand: 

&! = Q 1+2a; i / d 2 rr 1+2:r Vi(r) (40) 
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(it does not depend on Q) and 

B 1 = K 1 - 2 ^j J^ q ^j\ K ,q- K /2) (41) 
(it does not depend on k). Numerical calculations give 

b { p = 0.3145 e 2 Z 2 , b (L) = 0M377e 2 Z 2 , B 1 = 17.93 (42) 

Thus in the high-Q limit the expression for a 2 fully factorizes in the projectile and target. 
Its dependence on Q and k turns out to be intermediate between the eikonal and asymptotic 
triple pomeron predictions. It vanishes at large Q as l/Q 1+2xi , faster than the single pomeron 
exchange and asymptotic triple pomeron (~ 1/Q) but not so fast as the eikonal prediction 
1/Q 2 - It is also singular at k — > 0, but the singularity is weaker than predicted by the 
asymptotic triple pomeron. 

The double scattering amplitude following from (39) is 

M = -,CinV,^^ lBl ^ {Q *^ (43) 
where with a logarithmic accuracy (see (18)) 

d 2 K 

and we have introduced the nuclear radius Ra to make u>2 dimensionless. Passing to the 
low-x EMC effect we have to compare it with the dominant single scattering term 

Ai=iC\ lll2 {<S)P{s,ti) (45) 

where, calculated at k exactly equal to zero, 72(0) is finite: 



*2 = R\ +2X1 I 7^k~ 1+2xi Fa(k) In 2 k (44) 



72(0) = 2a s VN 2 - 1(R/N) (46) 

The EMC effect is characterised by the EMC ratio Remc of the nuclear structure function to 
A times the nucleon one. The double scattering contribution to it is given by Remc = 1 — 
where Xa = —A2/A1. From our formulas we find 

XA = C(A ~ 1)W2 R- A (QR A ) 2 ^ (^W 2 6XP [Ws ) (47) 
where the numerical coefficient is 

c = 2 9 - 2 ^7r 3 / 2 a 2 ci +) J B 1 (6 1T + b 1L )/(b T + b L ) (48) 

^,From (47) we expect the EMC effect to go to zero as Q —> 00 as Q~ 2xi , that is, rather 
slowly, essentially slowlier than one would find from the eikonal picture. Its dependence on 
A is enhanced by a factor ~ (Ra/R) 1 ~ 2xi \u 2 {Ra/ R) as comparted to the eikonal prediction. 
The enhancement is not so strong as one might naively predict on the basis of the asymptotic 
triple pomeron picture. 
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3 Some numerical estimates 

The obtained formulas in principal allow to calculate the double scattering contribution to the 
EMC effect on nuclei at small x. Since the EMC effect is small experimentally, one may hope 
that this contribution practically exhausts it (also see Conclusions for some justification). 
Our formulas contain two parameters: the strong coupling constant a s and the nucleon 
radius R. The value of a s can be extracted from the observed intercept A. As to R, it 
has to be of the order of the proton electromagnetic radius. Of course, one has to take 
into account inevitable uncertainties associated with the logarithmic character of the hard 
pomeron model and absence of scale in logs and logK. However in trying to apply our 
formulas to the experimental situation one meets with another serious obstacle. 

Comparison of the experimental proton structure function with the predictions based on 
the hard pomeron model indicates that the model can only be valid at very small x < xq = 
0.01. A power growth characteristic for the hard pomeron can be supported experimentally 
at x ~ x\ = 10~ 4 - 10~ 5 . The triple pomeron interaction picture requires that both s\ for the 
upper pomeron and s/s\ for the two lower ones be correspondingly large. It follows that this 
picture can only be applicable at extraordinary small x < xqX\ = 10 -6 . Besides,the large Q 2 
limit used to obtain a factorizable form for the contribution requires that {QRa) 1 ~ 2xi » 1. 
Present experimental data on the EMC effect do not satisfy these requirements. They are 
restricted to x > 0.005 and values of Q 2 below 2 (GeV/c) 2 for x < 0.01. Thus, comparison 
of our predictions to the existing data cannot be justified. 

Just for curiosity, if, notwithstanding this objection, one takes realistic values A = 0.2 
and R = 0.44 fm [9], then one gets for Ca (A = 40) at x = 0.0085, Q 2 = \A{GeV/c) 2 
A40 = 0.14 in an agreement with the experimental value A|§ = 0.154 ± 0.014. However with 
the same parameters one gets for Xe (A = 131) at x = 0.0065, Q 2 = 1.34 (GeV/c) 2 and C 
(A = 12) at x = 0.0055, Q 2 = 1.1 (GeV/c) 2 ) values A131 = 0.38 and \ 12 = 0.043 compared 
to the experimental values Affj = 0.16 ± 0.11 and Aff = 0.096 ± 0.013. 

Forgetting about the experimental situation, our predictions for the range of x and Q 2 
where one can expect the obtained formulas to be applicable are illustrated in the Table. 
One observes from it that the dependence on x is quite weak. In fact one obtains practically 
the same values of A in the whole range 10 -6 < x < 10~ 2 due to compensation of the growth 
with 1/x of both the numerator and denominator in (47). The predicted dependence on Q 2 
is clearly visible, although also not very strong as follows from (47). 
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4 Conclusions 

Two main consequences follow for the double scattering contribution to the structure func- 
tions of the nuclei from the hard pomeron model. First, at large Q 2 it behaves as Q~ l ~ 2xi 
where x\ = 0.3169 is a number which does not depend on the coupling and so is exactly 
known. As a result the corresponding low-x EMC effect should go to zero as Q~ 2xi . This 
property does not seem to be changed after the inclusion of higher order rescatterings. In- 
deed in the dipole approach one finds that they are described by the pomeronic fan dia- 
grams [4]. Then the initial pomeron coupled to the projectile will always be in the same 
regime as in our double scattering case with the only difference that the rest part will be- 
have as 5 nA for n rescatterings. The power scj governing the high Q 2 behaviour will be 
the smallest positive root of the equation (32) with 2A — * nA. This power grows with n: 
Xi = 0.3793, Xi = 0.4097, x^°^ = 0.5. Contributions from higher order rescatter- 
ings will then go down at high Q 2 faster than the double scattering contribution, although 
by a rather small power of Q 2 . Thus the behaviour Q~ 2xi = Q-0-6338 

seems to be a clear 

prediction for the small- x EMC effect which may serve to test the hard pomeron model. 

Second, due to singularities at small momentum transfers, the double scattering turns 
out to be more strongly dependent on A as compared to the standard eikonal predictions. 
This second prediction does not, however, seem so trustworthy, since it relies on the specific 
property of the lowest order hard pomeron model of having no intrinsic scale. 
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Table 

Predictions for the EMC ratios at largel/x and Q 2 . 



A 


X 


Q 2 (GeV 2 ) 


1 — Remc 


40 


10" 5 


10. 


0.0685 






100. 


0.0330 




10~ 6 


10. 


0.0775 






100. 


0.0374 


64 


10" 5 


10. 


0.104 






100. 


0.0502 




10~ 6 


10. 


0.118 






100. 


0.0568 


131 


10" 5 


10. 


0.190 






100. 


0.0914 




10~ 6 


10. 


0.215 






100. 


0.103 



